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ABSTRACT 


The  research  described  in  this  report  is  aimed  at 
characterizing  the  non-equilibrium  flow  and  coupled  chemistry 
associated  with  the  interaction  of  rocket  exhaust  effluents 
in  the  upper  atmosphere,  infrared  chemiluminescence  emitted 
from  the  reaction  between  0  and  C2H4  has  been  quantitatively 
measured  under  lean  conditions.  The  ratio  of  the  rate 
coefficients  for  the  reactions 

CHO  +  O  -»  CO  +  OH* 

and  CHO  +  02  CO  +  H02 

is  found  to  be  1.7  while  the  rate  coefficient  for  the  reaction 

CH20  +  0  -  CHO  +  OH 
-14 

is  found  to  be  less  than  3  x  10  cc/sec  with  a  quantum  yield 
of  less  than  20%. 

A  compilation  of  all  available  cross  section  information 
is  included  herein  for  use  in  highly  rarefied  plume  calculations. 

The  molecular  beam  measurements  on  the  interaction  potential 
between  O  atoms  and  plume  species  is  described. 
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PREFACE 


The  research  described  in  this  report  falls  into  four 
main  categories,  identified  in  the  contractual  work 
statement  as  follows: 

a)  Calculate  the  continuum  flow  structure  and  rate  chemistry 
of  high  altitude  plumes. 

b)  Calculate  the  vibrational  excitation  cross-sections  for 
plume  relevant  species  using  available  interaction  potential 
data. 

c)  Measure  the  quantum  yields  in  the  infrared  region  for  the 
chemi-excitation  involving  0  or  OH  and  plume  species. 

d)  Measure  the  inter-molecular  potential  between  atomic  oxygen 
and  various  plume  species  as  an  input  to  the  calculation 

of  vibrational  excitation  cross  sections. 

Item  (a)  is  discussed  in  Section  I,  item  (b)  in  Section  II, 
item  (c)  in  Section  III,  and  item  (d)  in  Section  VI. 
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I  -  Fluid  Mechanical  Studies 

In  this  portion  of  the  Wayne  State  effort,  emphasis  has 
beer,  focused  on  two  specific  problem  areas:  development  of  a 
1— D  unsteady  continuum  fluid  description  as  a  bridge  between 
the  2-D  steady  continuum  MULTITUBE  code  and  a  full  Monte  Carlo 
rarefied  flow  description  (under  development  at  TRW)  ,  and 
examination  of  the  methodology  for  including  chemistry  and 
excitation  processes  in  a  1-D  Monte  Carlo  rarefied  flow 
description. 

The  development  of  the  1-D  unsteady  continuum  flow  code 

has  been  described  in  the  previous  semi-annual  technical  report 

dated  15  December,  1971. ^  To  date  the  coding  on  this  program 

has  been  completed  and  the  calculation  successfully  checked  out 

for  the  case  of  expansion  into  a  vacuum.  Current  effort  is 

focused  on  studying  shock  formation,  i.e.  studying  the 

development  of  shock  structure  through  a  thick  shoe ^  description. 

This  phase  of  the  fluid  study  will  be  carried  on  by  7.  Boynton 
(2) 

elsewhere. 

The  Monte  Carlo  code  development  has  several  purposes  in 
addition  to  accounting  for  chemistry  and  excitatior  processes. 

These  include  an  examination  of  multicomponent  effects  in  a 
simple  geometry  and  as  a  means  of  gaining  familiarity  with 
the  structure  of  a  Monte  Carlo  analysis  for  advisory 
purposes  in  chemistry  input.  The  early  development  of  the  Monte 
Carlo  code  is  described  in  the  semi-annual  report  dated  15  December, 
1971,  so  will  not  be  further  described  here.  Work  in  this  last 
6  months  has  established  the  validity  of  the  1-D  unsteady  Monte 
Carlo  code  for  shock  wave  formation.  A  routine  for  the  inclusion 
of  vibrational  excitation  processes  has  been  developed  and 
included  into  the  M-C  analysis.  In  addition,  this  excitation 
subroutine  has  been  forwarded  to  TRW  for  inclusion  with  the  full 
Monte  Carlo  development  underway  at  TRW.  The  current  problem  in 
the  M-C  code  development  is  that  when  the  excitation  probability 
is  small,  the  code  is  slow  in  running  on  the  IBM  360/67  and  CDC 
7600  machines. 
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II.  Vibrational  Cross  Sections 


In  the  high  altitude  rarefied  flow  regime,  highly  non- thermal 
collisions  characterize  the  interaction  between  ambient  atmospheric 
species  and  exhaust  constituents.  In  order  to  characterize  the 
excitation  processes  associated  with  these  high  relative  energy 
collisions,  cross  section  information  for  various  plume  and 
atmospheric  species  must  be  generated.  The  close  coupling  code 
of  Marriott ^  has  been  applied  to  some  of  these  cross  section 
calculations  over  the  past  years.  in  this  reporting  period,  the 
results  of  Marriott's  past  cross  section  calculations,  in  addition 
to  some  new  results  of  plume  interest,  have  been  correlated  into 
a  form  convenient  for  evaluation  of  the  excitation  cross  section 
at  any  collision  energy  above  threshold. 

No  experimental  information  is  available  on  the  excitation 
cross  sections  for  any  of  the  collision  processes  calculated. 
Therefore,  in  order  to  attempt  to  assess  the  validity  of  these 
cross  section  calculations,  rate  coefficients  (or  relaxation  times) 
have  been  generated  by  integrating  the  cross  sections  over  a 
Boltzmann  distribution  of  relative  energies.  Where  thermalized 
data  are  available  these  calculated  relaxation  times  have  been 
compared  and  conclusions  drawn. 


Table  II a 


State  Notation  for  the  Polyatomic  Species 


Notation 


J 


A  full  description  of  the  theoretical  close  coupling  code  is 
currently  being  included  in  a  technical  report  describing  this  work^ 
and  is  also  available  elsewhere  so  only  the  results  will  be  summarized 
here.  The  species  considered  in  this  tabulation  include  CO,  C02 ,  H20, 
H2  and  He.  The  state  notation  for  the  polyatomic  species  is  given 
in  Table  n  a.  it  should  be  noted  that  the  v2  bending  mode  of  CO2 
is  doubly  degenerate  and  the  cross  section-  do  not  account  for  this, 
i.e.  the  total  excitation  cross  section  from  ground  to  the  vo  "  -*• 
state  is  twice  the  cross  section  listed  (see  Marriott  for  a  further 
discussion  of  this) . 

The  functional  form  assumed  for  correlating  the  energy  dependence 
of  the  cross  section  is 

"(2)  =  C3fc(E/Et^C2(ET/E)C4  n-i 


where  E  =  relative  collision  energy  (eV) 

E„,  =  threshold  energy 

1  2 
and  a ( E)  =  excitation  cross  section  (cm  ) . 

Table  lib  lists  all  the  coefficients  to  equation  (II-l)  for  the 
reactions  studied. 

These  cross  section  functions  can  be  integrated  over  a  Boltzmann 

distribution  of  relative  energies,  yielding  rate  coefficients  with  the 

inclusion  of  degeneracy  factors.  The  complete  tabulation  of  rate 

coefficients  as  determined  from  the  data  in  Table  lib  will  be  available 

in  a  forthcoming  report  complete  with  the  comparison  to  available  rate 

(2) 

coefficient  data.  The  conclusions  drawn  from  these  comparisons 

are  ,  that  at  high  temperatures  the  agreement  between  the  calculated  rate 
coefficients  and  measurements  is  good  in  all  cases  for  which  data 
exist  except  C02(v2)  +  H2 .  At  lower  temperatures,  i.e.  below  1000°K, 
the  calculated  rate  coefficients  are  generally  smaller  than  the  data 
suggest  reflecting  the  neglect  of  vibration-rotation  (VR)  and 
vibration-vibration  (W)  processes  in  the  calculation  of  the  cross 
sections.  These  latter  processes,  particularly  VR  processes,  are  expected 
to  be  important  in  effecting  thermalized  energy  transfer  processes 
involving  H2  and  H20.  For  the  case  of  C02(v2)  *-  H2>  the  calculated 


Table  II 


Coefficients  to  Cross  Section  vs.  Energy  Function 


REACTION 


N2  +  AR  Q(01) 

CO  +  HE  0(01) 

CO  +  HE  Q(02) 

CO  +  HE  Q(12) 

CO  +  H20  Q(01) 

CO  +  H20  Q(02) 

CO  +  H20  0(12) 

H20  +  H2  Q(01) 

H20  +  H2  Q(02) 

H20  +  H2  0(03) 

H20  +  H2  Q(04) 
H20  +  H2  0(12) 
H20  +  H2  0(13) 
H20  +  H2  Q(14) 
H20  +  H2  0(23) 
H20  +  H2  0(24) 
H20  +  H2  0(34) 
H20  +  H20  0(01) 
H20  +  H20  0(02) 
H20  +  H2C  0(03) 
H20  +  H20  Q(04) 
H20  +  H20  0(12) 
H20  +  H20  0(13) 
H20  +  H20  Q(14) 
H20  +  H20  Q{23) 
H20  +  H20  Q(24) 
H20  +  H20  Q(34) 
H20  +  C02  Q(01) 
H20  +  C02  Q(02) 
H20  +  C02  Q(03) 
H20  +  C02  Q(04) 
H20  +  C02  0(12) 
H20  +  C02  Q(13) 
H20  +  C02  Q(14) 
H20  +  C02  0(23) 
H20  +  C02  Q(24) 
H20  +  CO 2  Q( 34) 
C02  +  H2  Q(01) 
C02  +  H2  0(02) 
CO 2  +  H2  Q(03) 
C02  +  H2Q(04) 


ET 


2.89E-01 
2 . 660E-01 
5 . 310E-01 
2 . 660E-01 
2 . 660E-01 
5 . 310E-01 
2 . 660E-01 
1 . 980E-01 
3 . 950E-01 
4 . 410E-01 
4 . 650E-01 
1.980E-01 
2 .430E-01 
2 .680E-0I 
4 . 560E-02 
7 . 000E-02 
2 . 440E-02 
1.980E-01 
3 . 950E-01 
4.410E-01 
4 . 650E-01 
1 . 980E-01 
2 . 430E-01 
2 . 680E-01 
4.560E-02 
7 . 000E-02 
2.440E-02 
1.980E-01 
3 . 950E-01 
4 . 410E-01 
4.650E-01 
1.980E-01 
2 . 430E-01 
2 . 680E-01 
4.560E-02 
7.000E-02 
2 .440E-02 
8.280E-02 
1.660E-01 
1.640E-01 
2.910E-01 


C2 


1.093E03. 

3.658E00 

3 . 659E00 

3.779E00 

1.306E01 

1.175E01 

6.830E00 

3.937E00 

4.387E00 

2 . 628E00 

3.488E00 

2.398E00 

4.588E00 

5.631E00 

3.992E00 

4.180E00 

9.229E00 

5 . 160E00 

3. 3S8F.OO 

3.260E00 

3.134E00 

3.682E01 

3.977E00 

3 . 728E00 

6.307E0O 

5 . 153E00 

1.157E00 

1.098E01 

8.663E00 

1.481E01 

5 . 116E00 

4 . 923E00 

2.119E01 

1.451E01 

8 . 291E00 

7.518E00 

4.117E00 

3.748E00 

2.133EOO 

4.145E00 

5.473E00 


C3 


1.110E-23 
3 .263E-20 
8 . 189E-23 
1 . 299E-19 
1 . 449E-18 
2 . 189E-20 
4  .082E-22 
3.755E-16 
1 . 630E-16 
2 . 205E-18 
1.183E-18 
1 . 939E-16 
5. 263E-17 
2 . 619E-17 
5 . 143E-18 
3 . 589E-18 
5 . 393E-18 
2 .754E-17 

.  5. 112E-20 
2 . 065E-20 
1 . 676E-18 
2 . 756E-17 
5 . 328E-18 
1 . 993E-16 
1.067E-16 
1.419E-16 
1.708E-19 
9. 932E-21 
1.103E-19 
1.689E-23 
3 . 386E-19 
1.793E-16 
9 . 282E-18 
2.392E-17 
5 . 091E-17 
1.289E-16 
4.420E-16 
1.142E-16 
9.133E-17 
2 . 600E-17 


C4 


0.0 

-2.147E00 
-4.857E00 
-8. 742E-01 
3. 183E-00 
1 . 764E-00 
-2 . 357E-00 
1 . 380E-00 
1.734E-00 
-3 . 184E-01 
4.102E-01 
7 . 175E-01 
1.374E-00 
1 . 999E-00 
3.068E-01 
5 . 583E-01 
2 . 230E-00 
1.287E-00 
-3 .844E-01 
1 . 400^-00 
1 . 100E-00 
-1 . 709E-00 
1 . 144E-00 
9.232E-01 
2 .017E-00 
1.835E-00 
2 . 017E-01 
1 . 569E-00 
6 . 524E-02 
4.090E-00 
3 . 680E-00 
3 . 393E-01 
8.457E-00 
5 . 192E-00 
1.948E-00 
2 . 180E-00 
9. 704E-01 
1 . 244E-00 
3 . 772E-01 
1 . 176E-00 
1 . 721E-00 
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Table  con't. 


REACTION 

ET 

C2 

C 3 

C4 

C02 

+ 

H2 

Q(12) 

8 . 280E-02 

2.267E00 

8 . 757E-16 

9 . 479E-01 

C02 

+ 

H2 

Q(13) 

8 . 090E-02 

1. 395E00 

2 . 359E-17 

4 . 906E-02 

C02 

H2 

Q(14) 

2 . 090E-01 

3.476E00 

5.288E-18 

4 . 379E-01 

C02 

+ 

H2 

Q(  32) 

2 . 000E-03 

5.107E01 

8.2  54E-38 

7 . 115E-00 

C02 

+ 

H2 

Q(.24) 

1 . 260E-01 

7.104E00 

9. 580E-18 

1 . 683E-00 

C02 

+ 

H2 

Q(  34) 

1 . 280E-01 

4.391E00 

5 . 831E-18 

3 . 820E-01 

C02 

+ 

H20 

0(01) 

8 . 280E-02 

4.693E00 

1 . 256E-16 

1 . 173E-00 

C02 

+ 

H20 

0(02) 

1 . 660E-01 

5.237E00 

2 . 408E-17 

1 . 217E-00 

C02 

+ 

H20 

0(03) 

1 . 640E-01 

4.624E00 

1 . 387E-18 

2.679E-01 

C02 

+ 

H20 

0(04) 

2 . 910E-01 

4.696E00 

1.870E-20 

-1 . 638E-00 

C02 

+ 

H20 

0(12) 

8. 280E-02 

6.947E00 

5 . 245E-17 

1 . 795E-00 

CO  2 

+ 

H20 

0(13) 

8.090E-02 

6.495E00 

8. 598E-18 

1 . 726E-00 

C02 

+ 

H20 

0(14) 

2.090E-01 

1.164E01 

7 . 792E-16 

5. 282E-00 

CO  2 

+ 

H20 

0(32) 

2.000E-03 

1.030E01 

6. 608E-21 

1 . 440E-00 

CO  2 

+ 

H20 

0(24) 

1 . 260E-01 

1.488E01 

1 . 483E-15 

6 . 078E-00 

C02 

+ 

H20 

0(34) 

1 . 280E-01 

4.601E00 

9 . 341E-19 

-4 . 290E-01 

C02 

+ 

C02 

0(01) 

8. 280E-02 

1.955E01 

4. 993E-19 

4.957E-00 

C02 

+ 

C02 

0(02) 

1 . 660E-01 

1 . 144E01 

3 . 059E-20 

1 . 324E-00 

C02 

-i- 

C02 

n  /  m  \ 

y,  V  w  —  / 

1 . S40E-C1 

1  630E01 

A  ooor' 

O  COOn  r\  A 

C02 

+ 

CO  2 

Q(04) 

2 . 910E-01 

1.926E01 

1 . 209E-17 

7 . 240E-00 

C02 

+ 

C02 

0(12) 

8. 280E-02 

2.804E0C 

4 . 183E-18 

0.0 

C02 

+ 

C02 

0(13) 

8 . 090E-02 

1.064E01 

2 . 919E-19 

2 . 411E-00 

CO  2 

+ 

C02 

0(14) 

2 . 090E-01 

2.067E01 

1 . 880E-15 

8. 969E-00 

CO  2 

+ 

C02 

0(32) 

2.000E-03 

1.048E00 

2 . 062E-17 

0.0 

CO  2 

+ 

CO  2 

0(24) 

1 . 260E-01 

1.968E01 

7 . 164E-16 

7 . 333E-00 

C02 

+ 

C02 

0(34) 

1 . 280E-01 

8.858E00 

6 . 543E-20 

2 . 766E-01 

8 
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results  are  larger  than  measurements  suggest  even  without  the  inclusion 
of  VP  processes.  Apparently,  this  can  only  be  interpreted  by  assuming 
that  the  interaction  potential  used  in  the  calculation  is  not  quite 
correct.  This  point  is  being  currently  investigated. 

In  summary,  the  cross  section  information  listed  in  Table  II  b  is 
the  only  source  of  high  energy  excitation  data  available  for  use  in  the 
rarified  plume  -  atmosphere  chemistry  description.  These  data  are 
expected  to  be  reliable  for  this  application  even  though  the 
thermalized  comparisons  show  some  disagreement  with  data.  This  is  due 
to  the  expectation  that  at  high  relative  collision  energies,  VR  and  W 
processes  will  be  relatively  unimportant  compared  to  direct  VT  excitation. 

Current  work  on  the  close  coupling  code  is  aimed  at  including 
structure  in  the  collision  partner  in  order  to  characterize  W 
processes  for  species  of  plume  interest. 
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III.  Flow  Discharge  System  Measurements 

Prior  to  quantitative  model  predictions  on  the  radiative 
properties  of  exhaust  products  in  the  upper  atomosphere, 
several  important  chemical  reaction  rate  coefficients  must  be 
determined,  in  particular,  our  studies  have  focussed  on 
simple  exothermic  oxidation  reactions  involving  0  and  OH 
attack  on  various  hydrocarbon  fuel  species  and  accompanying 
vibrational  relaxation  processes. 

The  interest  simple  exothermic  reactions  is  directed 
toward  the  release  of  the  available  energy  of  the  reaction 
into  radiative  modes  of  the  product  molecules.  Measurements 
in  our  laboratory  have  been  made  from  2500  A  to  25  g,  in  an 
attempt  to  determine  the  reaction  mechanism  and  available 
radiative  modes. 


In  order  to  determine  which  of  the  various  hydrocarbon 

reations  was  worth  studying  in  detail,  i.e.  determining  quantum 

yields  at  2.7  p.  ,  a  brief  review  of  hydrocarbon  combustion 

chemistry  was  performed.  A  simplified  flow  diagram  for  the 

decomposition  of  straight  chain  hydrocarbons  is  shown  in  Fig.  hi  a  . 

We  note  that  O  atom  attack  on  cnH2n+2  Produces  OH  in  a  nearly 

thermo  neutral  reaction  in  addition  to  the  C  H~  radical.  This 

radical  is  rapidly  attacked  by  0  atoms  to  produce  CnH2n  and  OH 

with  sufficient  exothermici ty  to  produce  vibrationally  excited 
*  * 

OH  .  H20  can  subsequently  be  produced  by  the  reaction 

OH  +  OH.  This  cycle  continues,  as  outlined  in  the  figure  until 
C2H^  is  formed.  This  then  decomposes  as  shown  in  Fig.  Vlb  upon  further 
O  atom  attack. 


OXYGFN  ATOM  ATTACK  ■>/*  £TH'/££6'£ 


Figure  II  fb.  Schematic  of  the  Reaction  Mechanism  Proposed 


The  attack  of  0  atoms  on  C2H4  has  been  selected  as  a 
prototype  reaction  in  the  degradation  of  straight  chain  hydro¬ 
carbons  by  0  atoms.  From  Figure  Illb,  the  primary  products 

are  CHO  +  CH,.  The  rapid  reaction  of  0  +  CHO  is  sufficiently 
J  * 

exothermic  to  produce  OH  .  At  low  concentrations  of  C_H., 

*  ^4 

OH  is  expected  to  be  the  dominant  2.1  ^  radiator.  This  is 
discussed  at  length  in  a  recent  thesis  performed  at  Wayne  State 
University  under  the  direction  of  R.  Kummler. ^  The  work 
performed  during  that  thesis  study  will  be  published  as  a 
technical  report  under  the  plume  program. 

To  summarize  the  results  of  that  study,  experimental 
infrared  measurements  at  2.7^  were  taken  for  the  reaction  system 
illustrated  in  Figure  Illb  at  low  concentrations  of  C2H4< 

Under  these  conditions,  as  shown  in  the  computer  simulation 

* 

of  Figure  IIIc,  the  OH  is  considerably  higher  in  concentration 

★ 

than  H20  (both  of  which  radiate  at  2.7^).  The  kinetic 
mechanism  in  the  computer  model  is  shown  in  Table  Ilia.  Under 
these  lean  conditions,  strong  2.7u.  emission  is  observed  from 
the  primary  radiation  step 

CHO  +  0  -»  CO  +  OH*. 

The  primary  CHO  loss  step  is  found  to  be  the  reaction 

CHO  +  02  CO  +  H02- 

The  ratio  of  rate  coefficients  for  these  two  competitive  CHO 
reactions  is  experimentally  found  to  be  1.7.^  In  addition, 
the  rate  coefficient  for  the  slightly  exothermic  reaction 

CH20  +  0  -  CHO  +  OH 
-14 

is  found  to  be  less  than  3  x  10  cc/sec  with  a  quantum  yield 
of  than  than  20%.^  A  full  report  of  this  study  is  being 
prepared . 
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KINFTICS  OF  THE  REACTIONS. 
TIME  =  3. 65499E-03  SEC. 


NO. 

REACTION 

CONSTANT 

BATE 

1.  C2HU 

♦  0 

-CH3 

♦CHO 

8.UE-13 

6.CE 

14 

3.  HCHO 

♦0 

=CHO 

♦  OH* 

3.0E-15 

1.3E 

12 

6.  CHO 

+0 

=C02 

+  H 

3.0E-14 

7.4E 

11 

7.  CHO 

*H 

=C0 

♦  H2 

4.0E-12 

1  .  1  E 

13 

9.  CHO 

*02 

=C0 

♦  H02 

1.0E-12 

6.  9E 

14 

11  .OH 

♦OH 

=  H20 

♦  0 

2.5E-  12 

1.  IE 

11 

13.CH3 

♦  02 

♦  M 

=CH302 

♦  f! 

8.0E-32 

1.6E 

12 

15.  H20* 

♦  H 

=  H20 

♦  n 

1.CE-12 

1 . 7E 

12 

17. HCHO 

♦0 

=CHC 

♦  OH 

6.0E-  14 

2.  6E 

13 

19.  OH 

♦  HCHO 

=  H20 

♦  CHC 

1 . 4  E-  1 1 

1.5E 

13 

21.  OH 

♦  H2 

=  H20 

♦  H 

6.6E-15 

5.2E 

07 

23. HO 2 

♦0 

=  0H 

+  02 

1.0E-  11 

6.  9E 

14 

26. OH 

♦  HO  2 

=  H  20* 

♦  02 

1.0E-12 

1.7E 

11 

27. CHO* 

♦  0 

=co 

♦OH* 

1.CE-11 

8.0E 

10 

29.0 

♦  H02 

-OH  * 

♦  02 

3.0E-  14 

2.  IE 

12 

31  .CH302 

♦  HV 

-CH302 

*  HV 

1.0E-20 

7, 3E-11 

NO. 

2.CH3 

♦  0 

REACT  ION 

=HCHO 

♦  H 

CONSTANT 

3.0E-11 

BATE 

6.3E 

14 

4.  CHO 

♦0 

=C0 

♦OH 

1.0E-13 

2.  5E 

12 

6.  OH 

♦  0 

=02 

♦H 

5.0E-11 

9.  IE 

14 

8. OH* 

+  HV 

=0H 

♦  HV 

3.3E  00 

4 . 4  E 

10 

10. OH* 

♦  0 

=0H 

♦  0 

5 . 0  E-  11 

5.7E 

13 

12  .OH 

♦  OH 

=  H  20* 

♦  0 

2.5E-13 

1  .IE 

10 

14.  H 

♦  02 

♦  N 

=  H02 

♦  H 

3.0E-32 

2  .4  E 

13 

16.  H20 ' 

♦  HV 

=H2  0 

♦  HV 

3.  5E  01 

1.7E 

09 

18. CHO 

♦  0 

=C0 

♦  OH* 

2.0E-  12 

4.  9  E 

13 

20.  OH 

♦  HCHO 

=  H20* 

♦CHO 

1. 4E-12 

1  .5E 

12 

22. OH 

♦  H02 

=H20 

*0  2 

1.0E-11 

1.7E 

12 

24. H 

♦  H02 

=0H  * 

♦  OH 

3.0E-  15 

2.  3E 

10 

26.C2H4 

♦  0 

=CH3 

♦CHO* 

1 .0E-16 

7.2E 

10 

28.  H 

♦  H02 

=0H 

+0H 

l.OE-11 

7.7E 

13 

30. CO 

♦OH 

=C  0  2 

♦H 

1.5E-13 

1.7E 

11 

32.H20 

♦  C02 

=  H20 

♦C  02 

1.0E-20 

4.  2  E 

00 

Figure  Ilia. 


Kinetic  reactions  and  rate  coefficients  at  300°k 
for  the  computer  (Simulation  of  Figure  mb)  . 
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Experimental 

All  studies  were  performed  in  the  discharge  flow  system 
shown  schematically  in  Figure  IIId.Molecular  oxygen  (American 
Cryogenics,  purity  99.6%)  was  discharged  in  an  Evenson  or 
Broida  Cavity  at  2450  MHZ  using  a  Raytheon  PGM  10x2  diathermy 
system.  The  discharge  was  pulsed  at  238  Hz  for  simultaneous 
infrared,  visible  and  ultraviolet  measurements.  The  discharged 
oxygen  was  passed  into  the  flow  tube  reactor  through  a 
teflon  tube  with  a  radially-multiperforated  exit  stream  which 
permits  excellent  mixing  with  a  second  gas  and  prevents 
discharge  region  radiation  from  reaching  the  observation 
ports.  The  stainless  steel  reactor  is  lined  with  a  teflon 
sleeve  which  largely  prevents  wall  recombination  of  the  atomic 
species.  The  inside  diameter  of  the  reactor  is  one  inch  and 
typical  gas  velocities  are  30  to  40  m/sec  through  a  Heraeus 
Engelhard  E  -  225  mechanical  pump.  Six  ports  permit  radial 
viewing  through  KBr  windows  using  an  RCA  C31025C  uncooled 
photomultiplier  on  a  Jarrell  Ash  scanning  quarter  meter 
monochrometer  with  a  Keithley  602  Electrometer  recorded  on 
a  Honeywell  Electronik  194  recorder.  Signal  to  noise  ratios 
were  such  that  no  damping  or  time  integration  of  signals 
was  necessary  in  the  u.v.  and  visible.  Measurements  in  the 
infrared  were  made  using  a  Santa  Barbara  uncooled  PbS  detector 
(to  3^,)  ,  Ge: Au  detector  (to  8^)  and  Ge:Cu  detector  (to  25p,)  . 

All  infrared  measurements  were  made  with  a  PAR  HR8  Phase 
Sensitive  Lock  in  Amplifier. 


The  hydrocarbons  employed  as  reactant  gases  were  methane, 
ethane,  propane,  ethylene,  propylene,  trans-butene,  acetylene 
and  formaldehyde.  Linde  research  grade  methane  and  ethane 
were  employed.  Union  Carbide  chemically  pure  propane  was 
used.  The  propylene  and  trans-butene  from  Cryogenic  Sales 
had  purities  of  95%  and  99%  respectively.  Several  runs  were 
made  with  "house"  methane  with  no  perceptibly  different  results. 


Figure  Hid.  Schematic  of  discharge  flow  tube  reactor 
and  detection  system.  Legend: 

M  =  Monochromator  -  Jarrell  Ash,  quarter  meter, 
scanning . 

P  =  Photomultiplier  -  RCA  C31025C  (uncooled) . 

E  =  Electrometer,  D.C.  Amplifier  -  Keithley  602. 

Pi  =  Picoampere  source  -  Keithley  602. 

R  =  Recorder  -  Honeywell  Electronik  194. 

Hi  =  High  voltage  C.C.  power  supply  -  Power 
Designs  Pacific  model  2K-10. 


Summary  of  Wayne  State  Measurements 

A  summary  review  of  all  the  flow  discharge  measurements 

to  date  is  presented  in  Table  Illb.  It  should  be  mentioned 

that  reaction  (1)  in  the  table  is  the  generation  reaction 

for  producing  OH,  and  therefore  the  reactions  initiated  with 

OH,  i.e.  reactions  (2),  (3),  and  (4),  have  IR  signals  in 

excess  of  reaction  (1) .  Some  cf  these  preliminary  survey 

measurements  are  being  further  quantified  in  a  way  similar 

to  the  work  on  0  +  C_H.. 

2  4 
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Table  mb.  Wayne  State  University;  Summary  of  Data  to  June  1972 


CHEMILUMINESCENT  FLOW  TUBE  MEASUREMENTS 


Reaction 

System 

Visible 

2  /7y, 

5-10,* 

1. 

H  + 

NO--  NO  +  OH 

Meinel  bands 

large 

large 

(OH  + 

z 

H2-  H20  +  H 

signal 

signal 

OH  + 

OH-»  H20  +  0) 

2. 

OH  + 

CO 

nothing 

nothing 

? 

3. 

OH  + 

CH.0 

? 

large 

medium 

z 

signal 

signal 

4. 

OH  + 

C.H. 

nothing 

small 

medium 

/  z 

signal 

signal 

5. 

0  + 

C.H. 

Meinel  bands 

large 

medium 

Z  4 

signal 

signal 

6. 

0  + 

C4H8 

Meinel  bands 

large 

signal 

■> 

7, 

0  + 

C3H6 

Meinel  bands 

? 

? 

8. 

0  + 

CH-0 

Meinel  bands 

large 

? 

Z 

(3500-4500 

signal 

continuum) 

9. 

0  + 

C2H2 

Meinel  bands 

small 

signal 

? 

10. 

0  + 

C2H6 

Meinel  bands 

11. 

0  + 

CnH2n+2 

Meinel  bands 

12. 

0  + 

nh3 

HNO* 

Comments 

maximum  found  in 
2.7|i  signal,  no 
maximum  in  5-10  p, 


signals  greater 
than  (1)  above 


maximum  found  in 
2 . 7,j,  signal 


maximum  found  in 
2 . 7jj,  signal 
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IV.  Molecular  Beam  Experiments 

The  calculation  or  measurement  of  the  excitation  of  exhaust 
species  by  atmospheric  0  atoms  is  recognized  as  an  important 
problem.  To  date,  direct  experimental  measurement  of  these  excitation 
cross  sections  has  not  proved  feasible.  In  this  study,  the  excitation 
cross  sections  are  calculated  using  the  close  coupled  code  and  the 
intermolecular  potential  acting  between  0  atoms  and  plume  species.  j 

This  potential  is  being  experimentally  measured  using  an  Amdur-type  | 

molecular  beam  apparatus.  j 

This  experiment  proceeds  by  measuring  the  attenuation  of  a  1 

fast  (50  -  5000eV)  O  atom  beam  in  passing  through  a  scattering  gas. 

The  elastic  scattering  cross  section  is  then  determined  from  \ 

^(I/IQ)  =  -nS£  (VII.  1) 

where  I  =  attenuated  beam  flux 

I  =  incident  beam  flux 
o 

n  =  scattering  gas  density 
t  =  path  length  in  scattering  gas 
and  S  =  incomplete  scattering  cross  section. 

Thus,  from  a  plot  of  £n(I/I0)  vs. n  for  a  constant  i  t  the  cross  section 
can  be  determined.  For  a  mixture  of  2  atomic  species  in  the 
incident  beam,  i.e.,  0(3P)  and  0(^D),  the  (l/lD)  vs.  n  curve  can 
be  resolved  into  two  linear  parts,  yielding  two  cross  sections. 

S  is  an  incomplete  cross  section  since  the  finite  width  of  the 
detector  and  the  atomic  beam  cause  part  of  the  scattered  particles 
to  reach  the  detector. 

The  molecular  beam  apparatus  used  in  this  experiment  is  shown 
in  Figure  IV  a.  Oxygen  atomic  ions  are  produced  in  the  duo-plasmatron 
source,  focused  and  accelerated  through  a  magnet  to  separate  out 
both  neutrals  and  molecular  ions.  The  oxygen  atomic  ion  beam  is 
then  collimated  and  energy  selected  prior  to  being  neutralized  by 
charge  exchange  with  water  vapor.  The  neutralized  oxygen  atomic 
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ELECTROSTATIC 
DEFLECTION  PLATES 


Figure  IVa.  Amdur-type  molecular  beam  apparatus  at  Wayne  State. 


beam  is  then  passed  through  the  scattering  chamber  and  the 
attenuation  measured. 

The  beam  apparatus,  donated  to  Wayne  State  by  NASA,  originally 
had  problems  with  the  beam  handling  system.  The  original  uni¬ 
potential  lens  was  replaced  by  a  series  of  lenses  and  two  sets 
of  deflection  plates  to  guide  the  beam  through  the  small  orifices 
in  the  scattering  chamber.  The  system  now  is  operational  and 
preliminary  measurements  should  be  available  within  the  month. 

Measured  incomplete  scattering  cross  sections  depend  critically 
on  the  beam  geometry.  This  fact,  not  initially  recognized  in  the 
MIT  experiments,  is  displayed  in  Figure  IV  b.  Measurements  made 
too  close  to  the  cusp  in  the  scattering  cross  section  can  yield 
large  errors  to  the  sensitivity  in  the  cross  section  to  variations 

in  both  the  beam,  width  (d)  and  the  detector  width  (  d)  .  The 

.  2  5 

potential  determined  for  the  case  m  the  figure  is  V(R)  =  K/r  *  . 

The  repulsive  part  of  the  interaction  potential  determined 
from  these  beam  experiments  is  in  the  range  from  lev  to  about  30eV. 
No  difficulties  are  currently  anticipated  in  these  measurements. 
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Figure  ivb:  Calculated  dependence  of  relative  beam  width  ( d ) 
to  detector  width  ( p)  on  scattering  cross  section. 


V.  Rotational  Excitation 


The  rotational  excitation  of  simple  molecules  plays  a 
critical  role  in  determining  the  spectral  spread  of  near 
infrared  emission  from  vibrationally  excited  molecules.  Therefore, 
a  review  of  rotational  excitation  models  was  undertaken  to  provide 
background  in  estimating  the  spectral  signature  expected  from 
vibrational  emission  of  CO,  CC^  and  1^0  in  non-equilibrium  flow 
environments.  This  review  has  been  compiled  into  a  report  which 
is  in  the  final  stages  of  typing. ^ 

In  addition,  based  on  this  review  of  rotational  excitation 
processes,  short  numerical  codes  have  been  written  to  calculate 
the  effective  rotational  relaxation  times  of  simple  molecules, 
the  population  distribution  of  rotational  states  under  non¬ 
equilibrium  conditions  and  collision-limited  steady  state  emission 

from  rotational  states.  The  codes  are  based  on  Miller's 

(2) 

statistical  model  of  rotational  excitation  and  employ  the 
rigid  rotor  approximation. 
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